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A lumped-mass model has been used to represent the impact condition 
between a fuel duct and an IVHM (In- Vessel Fuel Handling Machine) . The non- 
linear effects of a Bellville spring and the free fall impact of the fuel 
duct on the IVHM were included. The purpose of the tests was to determine 
the loads on the fuel duct due to the impact. A comparison between experi- 
mental and theoretical re stilts is presented. 

INTRODUCTION 


In the operation of the FTR (Fast Test Reactor) an IVHM is used to handle 
core components under various component and core conditions. The component 
under consideration is a fuel duct which interfaces with the IVHM by means of 
a grapple mechanism. The grapple grips a nozzle -handling socket mounted on 
the end of the fuel duct. Under certain operating conditions the fuel duct 
drops a short distance to impact on the grapple fingers. A simulated mechani- 
cal mockup of the FTR core and the IVHM was set up, and IVHM fuel duct with- 
drawal and insertion tests were conducted. During these tests high accelera- 
tions were measured on the simulated fuel duct. The experimental data indi- 
cated that the acceleration loads were due to the impact of the fuel duct on 
the grapple. This paper describes how the results of the analysis and tests 
of the impact are interpreted in terms of the response of the simulated fuel 
duct and the IVHM structure. 

SYMBOLS 

a acceleration 

B scalar damper 

d drop distance 

g structural damping 

t time 


DESCRIPTION OF IVHM 


The FFTF (Fast, Flux Test Facility) simulated core test facility and an 
IVHM positioned for fuel duct insertion are shown in figure 1. The IVHM con- 
sists of a motor-driven mechanism mounted on a lead screw and a grapple 
assembly as shown in figures 1 to 3 - Four fuel duct support fingers (fig. 3 ) 
support the weight of the fuel duct during insertion and provide for the with- 
drawal forces. The IVHM is used to transfer fuel duct components from the 
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active core to fixed storage positions located outside the core region. Fuel 
handling problems can occur due to distortion of the individual fuel ducts and 
the core. The distortion effects can he induced by stainless steel swelling 
and thermal creep (ref. 1). Twelve core restraint yokes located hexagonally at 
two levels maintain the lateral position and compactness of the core. 


Free Fall Condition 

There are two operating conditions of the IVHM that have been investigated 
in this report. The first is a free fall impact of the fuel duct on the grapple 
and, secondly, a free fall impact combined with an initial condition of pre-load 
on the IVHM and fuel duct. During the insertion of the fuel duct, the above 
core load pad contacts the core (fig. 1), and a gap exists between the fingers 
on the grapple and the fuel duct. The core restraint yokes are retracted, the 
respective fuel duct position opens and the fuel duct drops the gap distance to 
impact on the grapple support fingers. 


Pre-Load Conditions 

Certain core conditions (ref. 1) require that a compressive force be 
applied to the fuel duct for insertion. The lead screw on the IVHM is used to 
apply the compressive force through the grapple (turnbuckle effect) to the fuel 
duct. Again, as the grapple engages the top of the fuel duct for the compressive 
load a gap exists between the grapple fingers and the nozzle -handling socket on 
the fuel duct. 

When insertion occurs the fuel duct experiences a free fall with an initial 
condition of pre-load. The compressive pre-load (turnbuckle effect) is relieved 
on both the fuel duct and the IVHM during the first milliseconds of fuel duct 
free fall. The fuel duct then impacts on the grapple support fingers. The load 
cell oscillogram indicated that the 1000 lb pre-load initial condition was a 
ramp function with a time duration of 0.026 second. The effect of the pre-load 
on the initial acceleration of the fuel duct is •unknown. Due to the manner in 
which the accelerations were recorded, the change in free fall time due to the 
pre-load could not be determined. 


LUMPED-MASS MODEL 


Basie Model 

After a review of the test facility and test data, it was decided to use a 
relatively simple model to represent the IVHM and fuel duct so that the results 
of the analysis would be available quickly. The physical properties used in the 
model in most cases were calculated from the results of load deflection and 
other tests on the simulated IVHM. Figure ^ shows the scalar elements that were 
chosen to represent the stiffnesses, masses, and damping properties of the IVHM 
and the fuel duct. The IVHM and fuel duct were idealized with six degrees of 
freedom with instrumentation locations and areas of nonlinear response considered 
in the selection of grid points and limped masses. The scalar damper values 
shown in the NASTRAN Input (fig. 5) were calculated from the frequency and 
measured structural damping values (see ref. 2 for equation). A structural 
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damping value (g = 15 percent) was calculated for the assembled IVHM by con- 
sidering the rate of decay and of successive rebounds of the fuel duct on 
grapple support fingers, as shown by the IVHM load cell oscillogram (fig. 6). A 
structural damping value of g = 2 percent was calculated for fuel duct. The 
fuel duct was empty except for a lead weight mounted in bottom of the duct to 
simulate the components of the fuel duct. 


Nonlinear Considerations 

Two nonlinear conditions were considered in the model shown in figure a 
Bellville spring with stiffness different in tension than in compression, 

Element k, and the free fall gap. Element 5* For these conditions the nonlinear 
provision of NASTRAN was used. The nonlinear effects were treated through the 
use of an additional applied load vector and by means of a transfer function, 
the NONLIN 1 and TF input cards. (See figure 5 for NASTRAN input.) To use the 
nonlinear provision requires rigid format 9, Direct Transient Response. In this 
format (ref. 3) numerical integration of the coupled equations is achieved by 
employing a central difference equivalent step-by-step procedure. It was found 
that the solution became unstable as the time step was increased or as the value 
of the lumped masses was reduced. The stability problem was evident through 
large and/or plus or minus oscillating displacements, A discussion of the sta- 
bility problems described here is found in appendix C of reference 4. 


COMPARISON OF RESULTS 


Results from the lumped-mass model for displacement, force, and acceleration 

5 3 

parameters are shown in figure 7 "to figure 12 for the 1-^, in. and in. free 

16 8 

fall conditions. On the two displacement curves (figs. 7 and 10) the theoretical 

free fall times ( d = — at^ ] for the 1-5- in. and 1— in. gap distances are shown. 

\ 2 J 16 8 


The theoretical time for a single mass provides a close approximation of the free 
fall time of the two degree of freedom fuel duct model. The peak displacement 
and acceleration did not occur at first time of impact. This is shown by both 
the NASTRAN and experimental results and is due to the response of the multi- 
degree of freedom system (figs. 6, J, and 10). 


A comparison of the wave forms between the experimental and NASTRAN results 
for the IVHM load cell (figs. 6 and 8) indicates that structural damping is 
larger than the 15 percent value used in the lumped-mass model. The structural 
damping (g = 15 percent) was determined from the rate of decay of successive 
rebounds. It was difficult to determine a consistent damping parameter after the 
impact event from the load cell oscillogram. 


It was found from the NASTRAN results that structural damping had a very 
small effect on the amplitude of the first impact event. However, the value of 
structural damping made a significant change in the response of the IVHM subse- 
quent to the first impact. Structural damping was varied from 0 to 15 percent. 
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The response of the IVHM due to the initial condition of the 1000 lb pre- 
load is shown in figure 11. The pre-load was input to the model at Grid Point 4 
as a ramp function with a time duration of 0.026 second. Wo pre-load was 
applied to the fuel duct Grid Point 5 for the reasons previously noted in this 
report. Nominally, the free fall time could he expected to be less than the 
theoretical time because of the increase in acceleration of the fuel duct. (See 
fig. 10. ) Nevertheless the response characteristics shown in figures 9, 10, and 
12 yield an approximation of the effect of the pre-load condition on the IVHM. 

Figures 9 and 12 show the acceleration response of the fuel duct and the 
IVHM. The bounce effect, due to the multi-degree of freedom response, and the 

successive rebounds are illustrated. In both instances, the 1^ in. and 1^ in. 

free falls, the accelerations on the IVHM are higher than on the fuel duct. 

A comparison (fig. 13) is made of the experimental and NASTRAN results for 
the accelerometer on the top of the fuel duct and the IVHM load cell. The 
effect of the pre-load is to reduce the amplitude of acceleration on the fuel 
duct but to increase it on the IVHM. (See figs. 9 and 12.) 

CONCLUDING REMARKS 


The experience gained in the application of NASTRAN to a complex reactor 
component impact condition has been presented. 

Comparison of the results obtained with NASTRAN with results from the 
experiments and- theory shows NASTRAN to be very effective in calculating the 
transient response of the IVHM and fuel duct. 

The fuel duct loads determined by these tests and analyses should be con- 
sidered qualitative. The simulated fuel duct has different dynamic properties 
than the prototype FTR fuel duct. The two mass models of the fuel drct should 
be extended to include the principal dynamic characteristics of the components 
of the FTR fuel duct. 
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FIGURE 3. IVHM SIMULATED GRAPPLE 
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FIGURE 4. 1VHM LUMPED -MASS MODEL 
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FIGURE 5. NASTRAN COMPUTER INPUT DATA 


NASTRAN EXECUTIVE CONTROL DECK 72 

ID GRAPPLE, RUN 102 

APP DISPLACEMENT 

SOL 9,0 

TIME 15 

DIAG 2,3,5,6,9,13 
CEND 

IVHM GRAPPLE IMPACT 

NONLINEAR BELLVILLE SPRING - 2 PERCENT STRUCTURAL DAMPING 


CASE CONTROL DECK ECHO 

CARD 

COUNT 

1 TITLE - IVHM GRAPPLE IMPACT TEST 

2 SUBTITLE - NONLINEAR BELLVILLE SPRING - 2 PER CENT STRUCTURAL DAMPING 


3 

SET 1 

» 3,4,5 

4 

SET 2 

• 5,6 

5 

SUBCASE 1 


6 

DLOAD 

• 1 

7 

TSTEP 

■ 10 

8 

NONLINEAR 

- 100 

9 

TFL 

* 1000 

10 

OUTPUT 


11 

LINE * 38 


12 

MAXLINES 

= 40000 

13 

NLLOAD 

* 1 

14 

ACCELERATION 

* ALL 

15 

DISPLACEMENT 

* ALL 

16 

ELFORCE 

= ALL 

17 

OLOAD 

• 2 

18 

OUTPUT (XYPLOT) 


19 

PLOTTER CALCOMP, 

MODEL 765,205 DENSITY 556 BPI 

20 

XAXIS 

= YES 

21 

YAXIS 

* YES 

22 

XGRID LINES 

■ YES 

23 

YGRID LINES 

=■ YES 

24 

XTITLE 

■ TIME SECONDS 

25 

YTITLE 

* DISPLACEMENT IN. 

26 

TCURVE 

= * * * • GRID PTS. 4 AND 5 ' 

27 

XYPLOT 

DISP / 4(T1) , 6(T1) 

28 

TCURVE 

- * * * * GRID PTS. 4 AND 6 ' 

29 

XYPLOT 

DISP / .401) , 601) 

30 

YTITLE 

* LOAD LBS. 

31 

TCURVE 

- NONL1N FORCE GRID PTS. 4 

32 

XYPLOT 

NONLINEAR / 4 01), 5 01) 

33 

TCURVE 

= LOAD CELL FORCE ELEMENT 3 

34 

XYPLOT 

ELFORCE / 3(2) 

35 

YTITLE 

» ACCELERATION IN./ SEC. 2 

36 

TCURVE 

* ACCELERATION GRID PTS. 4 

37 

XYPLOT 

ACCE /4 01) , 601) 

38 

BEGIN BULK 
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FIGURE 5. (cont.) NASTRAN COMPUTER INPUT DATA 

IVHM GRAPPLE IMPACT TEST 

NONLINEAR BELLVILLE SPRING - 2 PERCENT STRUCTURAL DAMPING 72 PAGE 2 




SORTED 

BULK 

DATA 

ECHO 





CARD 

COUNT 

. 1 . 

. 2 . 

. 3 . 

. 4 

. . 5 

. . 6 . 

. 7 . 

. 8 

9 . , 10 

1* 

CDAMP4 

101 

15.92 

1 






2* 

CDAMP4 

102 

535.03 

2 






3* 

CDAMP4 

103 

398. 

3 






4* 

CDAMP4 

104 

4.5 

4 






5* 

CDAMP4 

106 

44.26 

5 

6 





6* 

CELAS4 

1 

.040 +6 

0 

1 





7* 

CELAS4 

2 

2. 559*6 

1 

2 





8* 

CELAS4 

3 

. 714 +6 

2 

3 





9* 

CELAS4 

6 

.421+6 

*5 

6 





10* 

CMASS4 

11 

3.909 

1 






11* 

CMASS4 

22 

.453 

2 






12* 

CMASS4 

33 

.8022 

4 






13* 

CMASS4 

44 

.0854 

5 






14* 

CMASS4 

55 

.8420 

6 






15* 

CMASS4 

222 

.264 

3 






16* 

DAREA 

10 

5 


+1.0 





17* 

DAREA 

20 

6 


+L0 





18* 

DLOAD 

1 

1.0 

1.0 

2 

1.0 

3 



19* 

GRID 

100 


1.0 

.0 

.0 


23456 


20* 

GRID 

200 


2.0 

.0 

.0 


23456 


21* 

NOLIN1 

100 

3 


+L0 

200 

1 

200 


22* 

NOLIN1 

100 

4 


-L0 

200 

1 

200 


23* 

NOLIN1 

100 

4 


+1.0 

100 

1 

100 


24* 

NOLIN1 

100 

5 


-L0 

100 

1 

100 


25* 

TABLED1 

10 







TAB10 

26* 

+AB10 

ao 

ao 

0. 0005 

33.0 

10.0 

33.0 

ENDT 


27* 

TABLED1 

20 







TAB20 

28* 

+AB20 

0.0 

0.0 

0.0005 

336.0 

10.0 

336.0 

ENDT 


29* 

TABLED1 

100 







TAB30 

30* 

+AB30 

0.0 

0.0 

1.3125 

ao 

11.3125 

34.2+6 

ENDT 


31* 

TABLED1 

200 







TAB40 

32* 

+AB40 

-0.001 

-. 100+3 

.0 

0.0 

10.0 

42.4 +4 

ENDT 


33* 

TF 

1000 

100 

1 

+L0 




TF101 

34* 

+F101 

4 


+1.0 





TF102 

35* 

+F102 

5 


-1.0 






36* 

TF 

1000 

200 

1 

+1.0 




TF1Q3 

37* 

+F103 

3 


+10 





TF104 

38* 

+F104 

4 


- 1.0 






39* 

TLOAD1 

2 

10 



10 




40* 

TLOAD1 

3 

20 



20 




41* 

TSTEP 

ENDDATA 

10 

3000 

.0001 

5 
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FIGURE 6. EXPERIMENTAL RESULTS - 1 5/16 FREE FALL - TEST NO. 19. 5. 3 



5 


00 


X 

o 


o 

< 


a. 

oo 


o 



0 0.06 0.12 0.18 0,24 0.30 

TIME, SECONDS 












ACCELERATION, IN. /SEC 



TIME, SECONDS 








DISPLACEMENT, INCHES 


FIGURE 10. GRAPPLE AND FUEL DUCT DISPLACEMENTS - 1 3/8 IN. FREE FALL - 1,000 LB PRELOAD 
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FIGURE 11. LOAD CELL FORCE - 1 3/8 IN. FREE FALL - 1,000 LB PRELOAD 
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FIGURE 12 . GRAPPLE AND FUEL DUCT ACCELERATIONS - 1 3/8 IN. FREE FALL - 1,000 LB PRELOAD 



FIGURE 13. COMPARISON OF PEAK VALUES FROM EXPERIMENTAL 

AND ANALYTICAL RESULTS 


CONDITION 

LOAD 

CELL 

LBS 

ELEAAENT 
NO. 3 
LBS 

ACCELEROMETER 
TOP OF DUCT 
g's 

GRID POINT 
NO. 5 
g's 

15/16 IN. FREE FALL* 

4,420 

4,300 

49 

37 

* jk 

13/8 IN. FREE FALL 

4,580 

4,400 

32 

29 

* TEST NO. 19.5.3 






**1.375 IN. FREE FALL WITH AN INITIAL CONDITION OF 1,000 LBS 
(PRE-LOAD) - TEST NO. 19.8.2 



